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THE CHEMICAL PROPERTIES OF LOW-REDSHIFT QSOS 
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ABSTRACT 

We investigate the chemical properties of low-z QSOs, using archival UV spectra obtained with the HST 
and IUE for a sample of 70 Palomar-Green QSOs at z < 0.5. By utilizing the flux ratio of UV emission lines 
(i.e., Nv/C IV, (Si IV+O IV])/C IV, and N v/He II) as metallicity indicators, we compare broad-line region 
(BLR) gas metallicity with AGN properties, i.e., black hole mass, luminosity, and Eddington ratio. We find 
that BLR metallicity correlates with Eddington ratio while the dependency on black hole mass is much weaker. 
Although these trends of low-z AGNs appear to be different from those of high-z QSOs, the difference between 
low-z and high-z samples is partly caused by the limited dynamical range of the samples. We find that metal 
enrichment at the center of galaxies is closely connected to the accretion activity of black holes and that the 
scatter of metallicity correlations with black hole mass increases over cosmic time. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: abundances - galaxies: evolution - quasars: 
emission lines 



1. INTRODUCTION 

Measuring chemical properties of galaxies and their red- 
shift evolution is a crucial step in understanding galaxy evo- 
lution since the metallicity of galaxies is closely related to the 
history of star-formation, gas inflow, and outflow. A num- 
ber of observational studies have been devoted to measure 
the metallicity of galaxies, revealing metallicity correlations 
with various galaxy properties. In the local universe, it has 
been shown that metallicity scales with galaxy luminosity and 
mass (the luminosity-metallicity relation and mass-metallicity 
relation, respectively) based on the metall icity measured 
from gas emission lines dTremonti et al.ll2004[ and references 
therein) or stellar absorption lines (e.g.. lGallazzi et aLll2005b 
iPanter et al. 2008$). These scaling relations indicate that metal 
enrichment is closely connected to the galaxy mass assembly. 

Recently the redshift evolution of these scaling relations in 
star-forming galaxies has been extensively investigated. At 
z < 3, many studies have suggested apparent metallicity evo- 
lution as a function of redshift. The evolutio n is significant 
espec i ally for low ste l lar-mass galaxies (e.g.. ISavaglio et alJ 
2005 1 : |Erb et al.l l2006t [Maiolino et al.l 120081: IMannucci et ail 
2009t lYabe et al.ll2012f) , th ough it could be due to observa- 
tional selection effects (see IMannucci et aTj|2010l) . At z > 3, 
exploring the metallicity of galaxies is extremely challenging, 
because typical galaxies at such high redshifts are very faint 
and the classical metallicity indicators in the rest-frame opti- 
cal s pectra shift out of the atmospheric win dows (but see also, 
e.g. jLaskar et^l2^TTtlNagao et al.ll2012b . 

To extend the metallicity measurement toward higher red- 
shifts, one possible approach is to focus on active galactic 
nuclei (AGNs). Thanks to their high luminosity (Lagn) and 
a various metallic emission lines in their rest-frame ultravi- 
olet spectra, it is possible to infer the metallicity of broad- 
line regions (BL Rs) even for QSOs at z ~ 6 — 7 with ground- 
based telescopes (Kurk et al.l2007l:|juarez et al.l2009t see also 
iMortlock et al.ll201 ll) . Previously a positive relation between 
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the metallicity of B LRs (Zrtr) and the redshif t of AG Ns has 
been reported (e.g., iHamann & Ferlandl [l99l 119931) : how- 
ever it turned out that the apparent relation was caused by 
a selection bias and that the correlation between AGN lumi- 
nosity and BLR metall icity (the Lagn ~ Zb lr relation) was 
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fundamental (see, e.g., Nag ao et alJ 1200 6b). Note that the 
luminosity-metallicity relation of AGNs has been also re- 
ported based on the emission lines in the narrow-line re- 
gion (NLR), which is much more extended than BLR and 
traces th e chemical properties in the spatial scale of AGN host 
galaxies dNagao et al . 2006a; Mats uoka et al.ll2009l) . Interest- 
ingly, the luminosity-metallicity relation of high-luminosity 
QSOs shows no st rong redshift evolu ti on in the redshift 
range of 2 < z < 6 dNagao et al.ll2006allbl: Huarez et al.ll2009l: 
iMatsuoka et al.ll2.01 ibh . implying that the chemical evolution 
at the center of host galaxies is mostly com pleted at a very 
high redshift (see also lMatsuoka et aTJl201 1 al) . 

Although the observed Lagn _ Zblr relation in AGNs and 
its redshift dependence are crucial to constrain evolution- 
ary scenarios of the supermassive black hole (BH), the host 
galaxy, and the interplay between these two (i.e., the galaxy- 
BH coevolution), there are two main drawbacks that should 
be resolved. First, while the Lagn _ -Zblr relation is well 
established, it s orig in is still controversial. For example, 
IWarner et aTJ ( 120041) reported that the metallicity of BLR 
showed a correlation with the mass of BH (Mbh), but no 
correlation with the Eddington ratio (L/L^&S)- Their re- 
sult was recently conf irmed by a larger sa mple of QSOs 
( Matsuoka et alj[2TjTTbb . On the other hand, [Shemmer et alJ 
( 20041) claimed that the observed Lagn _ Zblr relation was 
caused by the dep endence of Zrt r on t he Eddington ratio, not 
on Mbh (see also iDietrich et al.ll2009l) . At lower redshifts, it 
has been reported that narrow-line Seyfert 1 galaxies (NLSls, 
whose Eddington ratios are believed to be high; see, e.g., 
iBorosonl 120021 ; iGrupd 120041) show higher me t allicity than 
typical broad-line AGNs (e.g..lWills et alj[l999t iNagao et ail 
120021: IShemmer & Netzeri 120021) . The higher metallicity of 
NLSls in these studies is qualitatively consistent with the re- 
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suit reported by Sh emmer et alj (120041) . 

The other drawback comes from observational limitations. 
The AGN metallicity based on both BLRs and NLRs has 
been examined predominantly for AGNs at z > 2, since the 
AGN metallicity studies generally utilizes emission lines in 
the rest-frame ultraviolet spectra. Consequently, the obser- 
vational studies with ground-based telescopes are limited to 
AGNs at z > 2. This prevents us from studying the differ- 
ence in the chemical proper ties of BLRs betwee n high-z and 
low-z QSOs. For instance, iNagao et al.1 (l2006bl) studied the 
£agn - Zblr relation of QSOs at 2.0 < z < 4.5. The pos- 
sible physica l origin of this Lagn - Zbtr relation has been 
exa mined by IShemmer et aT] (120041) for 2.0 < z < 3.5, and 
bv iMatsuoka et all d201 lbl) for 2.3 < z < 3.0. Although 
IWarner et al.l ( 12004 ) investigate this issue for AGNs in a wide 
redshift range of < z < 5, they did not examine the red- 
shift dependence of BLR chemical properties. Although there 
are a few attempts to infer the N LR metallicity based on 
the rest-frame optical spectra (e.g.. IStorchi-Bergmann et al.l 
119981: INagao et al.ll2002l iGroves et al.ll2006l) . those methods 
are in turn difficult to apply for high-z AGNs since infrared 
spectroscopy is required to obtain rest-frame optical emis- 
sion lines, thus inconvenient for the comparative study be- 
tween high and low reds hifts. Using space observati ons with 
the IUE and the HST, IShemmer & Netzeii (12002b investi- 
gated UV spectra of low-z AGNs and showed that a signif- 
icant Lagn - Zblr relation was present also in low-z AGNs. 
However, they did not examine the physical origin of the 
£agn - Zblr relation and thus it is not clear whether BH mass 
or the accretion rate derives the observed Lagn _ Zblr rela- 
tion, and whether the low-z and high-z AGNs show the same 
correlations. 

Motivated by these considerations, in this paper we inves- 
tigate the correlation between BLR metallicity with various 
AGN properties, including mass, luminosity, and Eddington 
ratio, for a sample of 70 low-z Palomar-Green (PG) QSO at 
z < 0.5, by utilizing the archival UV spectra. We describe the 
sample selection and the data in §2, the data analysis and the 
fitting procedure in §3. The main results are presented in §4, 
followed by discussion in §5 and summary and conclusions in 
§6. We adopt a cosmology of Ho = 70 km s -1 Mpc -1 , 57a = 0.7 
and f2 m = 0.3. 

2. SAMPLE SELECTION AND THE DATA 

2.1. Sample selection 

To investigate metallicity of high-z QSOs, both permitted 
and weak semi-forbidden lines have been used. The weak 
semi-forbidden emission lines in the rest-frame UV spectra, 
i.e., N iv]A1486, O IH1A1 663, and N mlA1749 are good 
metallicity indicators (e.g.. IShieldslll976t [Baldwin & Netzeii 
1 1 978t lOsmenl 1 980t lUomotoll 1 984t IWarner et al.ll2002l) . since 
the flux ratio among these lines do not show strong depen- 
dences on physical properties of gas clouds (such as the den- 
sity and ionization parameter). However, the actual appli- 
cation to observational data is generally difficult since those 
semi-forbidden emission lines are too faint to be measured 
accurately. Thus, stronger emission lines are preferred in the 
studies of Zblr- For example, the flux ratios of N VA1240 
to C IVA1549 and N VA1240 to He IIA1640 have been uti- 
lized to infer Zblr by comparing them with the prediction 
of photoionization mode l s (e.g.. [Hamann & Ferlandl 119921 
1993 ] iFerTaiid et alj|1996t iKorista et al 1 1 19981: IDietrich et alj 
1999t IDietrich & Wilhelm-Erkensll2000l; Hamann et al.ll2002b 



IDietrich et aT] 120031) . INagao et al.l d2006bl) showed that the 
flux ratios of (Si IVA1397+0 iv]A1402)/C IVA1549 and 
Al IIIA1857/C IVA1549 are also useful to infer Z B lr through 
photoionization m odel runs, that are used for inferring Zbi ,r i n 
high-z QSOs (e.g. Jjuarez et al.ll2009l:lMatsuoka et alJ20l'lb|) . 

Since most previous studies used the rest-frame UV spec- 
tra to infer BLR metallicity of high-redshift Type-1 QSOs, 
we therefore focus on high-luminosity QSOs at low-redshfit, 
for which UV spectra are available, in order to investigate the 
chemical properties of low-z QS Os compared to high-z QS Os. 
The Palomar-Green (PG) QSOs (ISchmidt & GreerJI 19831) are 
well-studied low-z luminous Type-1 AGNs, and the UV spec- 
tra of many PG QSOs have been previously obtained with 
space facilities, thus suitable for our BLR metallicity study. 

We selected all PG QSOs at z < 0.5 (89 objects), for which 
reliable black hole masses are available from either rever- 



beration mapping results dPeterson et al.l 12004c iDennev et al 



2010) or single-epoch estimates (IVestergaard & Peterson 



2006). To investigate BLR metallicity, we will use the 



flux ratios of N VA1240/C IVA1549, N vA1240/He HA1640, 
and (Si IVA1397+0 IV]A1402)/C IVA1549, since these emis- 
sion lines have relatively large eguivalent widths. Thus, we 
searched for available UV spectra previously obtained from 
space facilities, using the Mikulski Archive for Space Tele- 
scopes (MAST). Among 89 PG QSOs at z < 0.5, the archival 
UV spectra covering the reguired emission lines were avail- 
able for 86 objects. Among them, we excluded 7 broad ab- 
sorption line (BAL) QSOs, since the strong absorption fea- 
tures prevent us from measuring the emission-line fluxes ac- 
curately. We also excluded 9 additional objects, for which 
the spectral guality is too low to identify the aforementioned 
emission lines. Thus, we finalized a sample of 70 PG QSOs 
for this work as listed in Table 1 . 

2.2. Data 

We obtained the UV spectra taken with International Ul- 
traviolet Explorer (IUE) or Hubble Space Telescope (HST) 
through the MAST database. We collected all available spec- 
tra of our targets and used the best guality spectrum with an 
order of Cosmic Origins Spectrograph (COS), Space Tele- 
scope Imaging Spectrograph (STIS), Faint Object Spectro- 
graph (FOS), and IUE when multiple instruments have been 
used. In summary, we utilized 10 COS spectra, 4 STIS spec- 
tra, 26 for FOS and 30 IUE spectra. 

Specifically, we used the SWP (1200 -2000 A) data of IUE, 
G130H (1150- 1600A) and G160H (1600 -2300 A) data of 
HST/FOS, and G140M (1150- 1740A) data of HST/STIS, 
for lower-z QSOs. For relatively higher-z QSOs, we used the 
LWP ( 1 800 - 3200 A) data of IUE and G270H (2300 - 3200A) 
data of HST/FOS data. Finally, for the H 577COS data, we 
used the combined data in two spectral ranges, i.e., G130M 
(1 150- 1450A) and G160M (1405 - 1775A), in order to cover 
the N V and C IV lines at the same time. 

We combined the spectra of each exposures by calculat- 
ing the error-weighted mean. For the STIS data, we com- 
bined the spectra using the exposure-time as a weight, be- 
cause we could not eliminate artificial spark features effec- 
tively when we adopted the error-weighted mean. However, 
the spectra are gualitatively consistent. In the case of COS 
spectr a, we used IDL rout ines developed by the COS GTO 
team ( iDanforth et alj|20lol) . We smoothed the spectra in the 
wavelength direction by adopting 7 pixel smoothing for the 
COS data and 2 pixel smoothing for the STIS data. 
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TABLE 1 

Log of archival UV data & AGN properties 



Object 


Redshift 


Observation Date 


Telescope/Instrument 


log[M BH /M ] 


Ref. 


log[L bo i/ergs M 


S/N 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


PG0003+158 


0.450 


1 nm m nc f\i 

1993 Nov 05,07 


Tin Ti 11 1 / ^\ £1 

HS 1/POS 


n i c i nno 

9.25 ±0.03 


2 


A f* £ A 1 n 1 o 

46.64 ±0.1 8 


5.47 


P(j0003+iyy 


0.026 


ani n l 1 1 no 

2010 reb 08 


HS 1/COS 


*7 1 o 1 n 1 1 

7. 13 ±0.11 


1 


^ a a o 1 n n /i 

44.43 ± 0.04 


28.09 


tip nnn "7 i 1 f\£. 

PG000/+106 


0.089 


1 no 1 t no 

198 1 Jun 0o 


TT TT7 /CU (n 

lUb/SWP 


o t 1 in aa 

8.71 ± 0.09 


2 


a a i n i n 

44.96 ± 0.10 


7.30 


PG0026+129 


0.142 


| flA A \[ 1 i-i 

1994 Nov 27 


Tin rp 11 1 / ^\ £1 

HS1/POS 


o ct 1 n 1 1 

8.57 ± 0.1 1 


1 


c ci i n no 

45.57 ± 0.03 


34.59 


PC0049+ 1/1 


0.064 


1985 Jul 3 1 


TT TT7 /C 

iur/a WP 


o 1 1 _l n no 
8.33 ± 0.09 


o 
2 


A A 1£. 1 a n 

44.36 ± 0.23 


3.22 


Drrift^n i 1 O/i 
PCj0050+124 


0.061 


19/9 Dec 22,23 


TT TTI7 IQWJT} 

lUb/aWr 


h /ii _i_ n i n 
/.42 ± 0. 10 


o 

2 


A a £.t\ _L n n^ 
44.69 ± 0.06 


12.12 


tip nnn ■ ic i 

PG005 2+251 


0. 155 


i Ann t in 

1992 Jun 29 


TT TT7 7Ti 


o cc i n nn 

8.55 ± 0.09 


1 


A C OO _l_ A nc 

45. /8 ± 0.05 


14.40 


PCjOI 5 /+UU1 


n l (*a 
0.1o4 


i noc a no 
1985 Aug 09 


TT TTJ /C\XfD 


q i c i n no 

8.15 ± 0.09 


o 

2 


a c on j_ n (yi 
45. /0 ± 0.0 / 


10.30 


PG0804+761 


0.100 


inin t 11 

2010 Jun 12 


TTPT /P C" 


o oo I n nc 

8.82 ±0.05 


1 


/i c nn I n ni 

45.99 ± 0.02 


62.65 


PG0838+7 /U 


0.131 


innn c i/i 

2009 Sep 24 




O 1 i in AA 

8.13 ± 0.09 


2 


A C 1 A _l_ A ft£ 

45.24 ± 0.06 


19.15 


PC0844+34y 


0.064 


198/ JNOV 30, Dec 01 


TT TTJ /C\XfTJ> 


t nc _i_ n i ° 
/.y5 ± 0. 18 


1 


A c nn _L n n*^ 
45.00 ± 0.06 


12.03 


np Afti 1 i cic 

PO0921+525 


0.035 


inoo tt^u io in 

1988 reb 28,29 


TT T~TT /CAT m 

lUh/SWP 


*7 1 O 1 A 1 1 

/.38 ± 0.1 1 


1 


a A in I n nc 

44.30 ± 0.05 


14.94 


dp nrn i i 1 in 

PO0923+129 


0.029 


inoc ~\ t . ■ n 1 

1985 May 01 


TT TT" /CAT fn 

lUhnSWP 


o c o I n 1 n 

8.58 ± 0.10 


2 


A A i a in no 

44.14 ± 0.08 


9.91 


dp An /i t i on/: 

PCj094 /+396 


0.206 


i nn/i a ii n/C 

1996 May 06 


T TOT /T7/~\C 

HS l/r(JS 


/ i in A A 

8.66 ± 0.09 


2 


/l C O A 1 A 1 /I 

45.84 ± 0.14 


6.78 


i")/" i i a i i r\Ar\ 

PCjIOI 1-040 


0.058 


inin Air i/; 

2010 Mar 26 


t tot /rrio 
HS 1/COS 


*7 i n i n nn 

/.30±0.09 


2 


aao') I n m 

44.83 ± 0.03 


35.03 


II/- ^ iaii , AAO 

PCj 101 2+008 


0.185 


i nnn a . in 

1990 Apr 10, 


TT TT" /CTT fl) 

lUb/SWP 


O 11 I A A A 

8.23 ± 0.09 


2 


A c A i I n 1 n 

45.41 ± 0.10 


7.46 


11/' i i nil . r i n 

PG 1022+5 19 


0.045 


1 AO") \yl TIT ni 

1983 May 31;Jun 01 


TT TT" /C\T / 11 

JUb/SWP 


f n 1 i n 1 n 

6.32±0.19 


1 


a A a o 1 n nn 

44.48 ± 0.09 


8.53 


PCj 1048+342 


0.167 


i nn t m 10 

1993 Nov 13 


TT TT" /CAT 7"A 

iUb/SWP 


o o c i n nn 

8.35 ±0.09 


2 


A A 1 A 1 n /" o 

44.74 ± 0.68 


1.12 


PCj 1049-005 


0.357 


1 nni a n i 

1992 Apr 01, 


T TCT /T7/~VC 

HS 1/rOS 


n i / in nn 

9.16±0.09 


2 


a s~ n i n oi 

46.17 ± 0.22 


4.36 


dp 1 1 m nn/: 

PCj 1 1 03-006 


0.425 


1 nni t~\ in 

1992 Dec 29 


T TCT" /T7/~VC 

HS 1/rOS 


A i a I n 1 n 

9.30±0.10 


2 


<A 1 1 1 A 1 1 

46.1 1 ± 0.1 1 


8.68 


n/^ i i i c . a r\i 

Pul 1 15+407 


0.154 


i nn/"" ii it 1 n 

1996 May 19 


T TCT" /I "/ "\ C 

HSl/rOS 


*i f c i n nn 

7.65 ±0.09 


2 


i c /n i nno 

45.62 ± 0.08 


12.30 


rip | 11 £ , t 1 c 

PCj 111 6+2 1 5 


0. 177 


i nni c-iU i n on 
1993 reb 19,20 


TJCT /T7p. C 

HS 1/rOS 


o c i _i_ n nn 
8.51 ± 0.09 


o 
2 


a in -j_ n i c 
4O.30 ± 0. 15 


a ca 
6.56 


HP i i in i 1 in 

PCjI 1 19+120 


0.049 


i noi m. o 1 izT 

1982 JNov 21,26 


TT TT~I /CAT TTi 

lUh/SWP 


n A c in nn 

/.45 ± 0.09 


2 


/I /l £1 _l_ A (\£L 

44.62 ± 0.06 


12.89 


np i m , a oo 
PCjI 121+422 


Pi ^1A 
0.234 


1 ftftC T n „ AO 

1995 Jan 08 


TT TC /QWTD 

lUb/bWr 


o ni _i_ n nn 
8.01 ± 0.09 


o 
2 


a c n/i _l n i i 
45.94 ± 0. 1 1 


7.06 


DP 1 | iin 1 1 A 

PCjI 14y— 1 1U 


n n/ln 
0.049 


i nno on 
1992 Dec 29 


TT TTJ /QWTD 


-7 nn i n i n 

/.yo ± 0. 10 


o 
2 


/i /i o c _l n i n 
44.25 ± 0. 10 


7.58 


DP 1 1 C 1 i 1 1 T 

PCjI 151+1 17 


0.176 


i ao7 t in in 

198 / Jan 29,30 


TT Tn /CAT TTi 

lUh/SWP 


o ci i n nn 

8.53 ± 0.09 


2 


<c /"c l A 1 1 

45.65 ±0.11 


6.68 


n/" 1 1 ini , aoi 

PCj 1 202+281 


0.165 


i nn/"" t ill 

1996 Jul 21 


T TCT" /I "/ "\ C 

HSI/rOS 


o ca i n nn 

8.59 ±0.09 


2 


A A f\C \ A 1 ^ 

44.95 ±0.14 


7.07 


DP 1111 i 1 A 1 

PCjl21 1+143 


0.085 


inni T7„u r\A m 

2002 reb 04,0/ 


T TCT" /CT"TC 
HS1/S11S 


*7 n^ i n a a 

/.94 ± 0.09 


2 


A C C 1 1 f\ f\A 

45.63 ± 0.04 


19.41 


dp i o 1 £ i n^n 
PC1216+069 


0.334 


i nni w in,- i c 
1993 Mar 15 


HS 1/rOS 


o i q A- n no 
9. 18 ± 0.09 


o 
2 


A& CO l n 1 A 

40.52 ± 0. 14 


7.03 


dp i ii£ i nn 

PCj 1226+023 


0. 158 


i nn i t. . i n 

1991 Jul 9 


T TCT /T7PC 

HSI/rOS 


o ni i n aa 

8.93 ± 0.09 


1 


a £. ca I a nn 

46.59 ± 0.09 


10.60 


1 1 /■ i inn i i n /I 

PCj 1229+204 


0.063 


i n o i \ a ni ni 

1982 May 01,02 


TT TT" /CAT / l"» 

IUb/SWP 


'H O A 1 nil 

7.84±0.21 


1 


a c 1 1 I n n /I 

45.11 ±0.04 


21.02 


dp 1 1 /i /i i ni /: 

PCj 1244+026 


0.048 


i ao T7„u no 

1983 reb 08 


TT Tn /CAT TTi 

IUb/SWP 


c cn I n nn 

6.50 ±0.09 


2 


h /i ">n I n 1 n 

44.30 ± 0.10 


7.72 


dp iicn i cm 

PU1259+593 


0.472 


i nn i oo 

1991 Dec 2/ 


T TCT /T7/~\ C 

HSI/rOS 


o nn I n i n 
8.90 ± 0.10 


2 


a f i f I n n 

46.76 ± 0.23 


4.15 


t^/^ 1 ~i ai 1 no 

PCjI 302-102 


0.286 


ino/"" t l ic i s 

1986 Jul 25,26 


TT TT" /CAT TTi T AT TTi 

IUb/SWP, bWP 


o o f 1 n 1 n 

8.86±0.10 


2 


a f a c 1 n n 1 

46.45 ± 0.01 


57.25 


n/" 1 1 1 n^7 i no c 

PG1307+085 


0.155 


1 AOA 11 f\ A 

1980 May 04 


TT TT" /CAT TTi 

IUb/SWP 


o /^i 1 n 1 1 

8.62±0.12 


1 


A c on 1 n n^7 

45.80 ± 0.07 


10.33 


1 1 /■ i i "i 1 n i no 

PO1310-108 


0.035 


i nn c t -1 i_ i i 

1995 reb 1 1 


TT TT" /CAA TTi 

IUb/SWP 


I I n nn 

7.86 ±0.09 


2 


a A 1 1 1 n 1 n 

44.13 ± 0.10 


7.66 


DP 1 111 i llCA 

PCj 1322+659 


0.168 


i nm t 1 a 

1 99 / J an 19 


T TCT /T7/~VC 

HS 1/rOS 


o iz; I n 1 1 

8.26 ± 0.1 1 


2 


45.52 ± 0.04 


23.02 


DP IT/II i 1 C O 

PCjI 341+258 


0.087 


i nn c a it 11 

1995 Mar 22 


TT TT~I /CAT •' IT 

lUh/SWP 


o ni I n 1 n 

8.02 ±0.10 


2 


A A 1 1 I A 1 1 

44.71 ± 0.13 


5.70 


dp nci i £ n c 

PCj 135 1+695 


0.030 


im i t n 

201 1 Jun 27 


T T C T /P C 

HS 1/COS 


i ci I nil 

7.52 ±0.12 


1 


/to zc i I n 1 n 

43.63 ± 0.10 


8.61 


DP 1 1 CI i 1 Ol 
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53.54 


PG1617+175 
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45.24 ±0.09 


8.66 


PG1626+554 
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8.48 ±0.09 
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8.98±0.10 
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46.25 ±0.18 
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2010 Oct 28, 
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0.89 
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NOTE. — Col. (1): Target ID. Col. (2): Red shift. Col. (3): Observ ed date. Col. (4 ): Telescope and Instrument. Col. (5): Black hole mass from [Peterson et al] 120041) : 
|\^stergaai'd & Peterson i 2006); Denney et al. 1 201(3) with a new virial factor ( Woo et al. 20Tc|). Col. (6): References for reds hift and black hole mass. 1 - R everberation -mapped AGNs 
i Peterson et al . 2004), 1* - Reverberation-mapped AGNs iDenney et aU20ld) , 2 - AGNs with single-epoch black hole mass fVestergaard & Peterson 2006). Col. (7): AGN bolometric 
luminosity calculated from the monochromatic luminosity at 1350A by multiplying a bolometric correction factor, 3.81. Col. (8): Signal-to-noise ratio per resolution element at 1350A 
in the rest-frame. 
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If a target has been observed at multiple-epochs, we chose 
only one epoch with the best data quality to avoid any time- 
variation effects. Table 1 lists the observation data and the in- 
strument for each target. In this table we also list the signal-to- 
noise ratio per resolution element calculated at the rest-frame 
1350A continuum. 

2.3. AGN properties 

To compare with BLR metallicity, we measure and col- 
lect other AGN properties, i.e., black hole mass, bolomet- 
ric luminosity, and Eddington ratio. We collected black hole 
mass of the sample QSOs, which has been previously deter- 
mined by the reverb e ration mapping stud ies for 18 objects 
dPeterson et al.ll2004t iDennev et ail 1201 01) or by the single - 
epoch method for 52 objects JVestergaard & Petersonll2006l) . 
A black hole mass measurement based on reverberation- 
mapp ing results is availble for PG121 1+143 dPeterson et al.l 
120041) . however it has large uncertainty due to the low data 
quality, and it has been excluded in other reverberation 
sample studies. Thus, we will use single-epoch mass for 
PG1211+143. 

We re-calculated black hole mass of the sample by adopt- 
ing the updated virial factor of 5.2 (IWoo et al.ll2010l) . which 
is slightly small er than the previous v irial factor (5.5; 
lOnken et al.ll2004l see also lPark et aD2012l) 

As the uncertainty of blac k holes mas ses, we 

adopted the values given bv iPeterson et all (120041) : 
IVestergaard & Petersonl (120061) ; lDennev et al.l d2010l) . 

For AGN bolometric luminosity (LagnX we use d the ob- 
tained UV spectra to measure monochromatic luminosity at 
1350A, which is presumably not heavily contaminated by the 
host galaxy stellar light. To measure the flux at 1350A, we 
fitted the AGN continuum between 1210Aand 1700 A with 
a power-law function. The measured 1350A monochromatic 
continuum luminosity is then used for calculating AGN bolo- 
metric lu minosity by multi plying a bolometric correction fac- 
tor, 3.81 dShen et al.ll2008l) . Note tha t this bolometric correc- 
tion factor is the same as adopted bv iMatsuoka et al.l (1201 lbl) 
for high-z QSOs. The measurement uncertainty of AGN lu- 
minosities was calculated based on the signal-to-noise ratio of 
the spectra. 

In Figure 1, we present the distribution of the sample prop- 
erties; the redshift, black hole mass, bolometric luminosity, 
and Eddington ratio. The black hole mass ranges over 3 or- 
ders of magnitude (from 6.32 to 9.71) with an average of 
8.26 ±0.71 Mq. The bolometric luminosity also ranges over a 
large range from 10 43 6 to 10 46,8 erg s -1 while some fraction of 
the sample has relatively low luminosity and belong to Seyfert 
class rather than QSOs. The mean Eddington ratio of the sam- 
ple is ^10% with 0.66 dex dispers ion, indicating that there is 
a large range of accretion activity dWoo & Urr v 2002). 

3. ANALYSIS 

3.1. Multi-component fitting procedure 

The flux measurement of BLR emission lines is an im- 
portant step for investigating Zblr- It is known that BLR 
emission lines some times show s i gnificant asymmetric ve- 
locity profiles (e.g.. |Corb inl [l997t IVa nden Berk et all l200lt 



iBaskin & Laorl 120051) . hence a single-Gaussian model does 
not generate a reliable fit for such cases. To fit asym 
metric velocity profiles of QSO UV emission lines, vari 
ous models have been adopted. Here, we examined 4 mod 
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FIG. 1 . — The distributions of the redshift, BH mass, AGN luminosity, and 
the Eddington ratio of the sample. 

els, namely, double-Gaussian, Gauss-Hermitian, modified- 
Lorentzian, and 2 power-law functions to determine the best 
line profile to use. In Figure 2, we present the C IV line 
of PG0003+158 as well as 4 different model fits, which 
show slightly different results, particularly at the wing of the 
line. Through our visual inspection, we decided to adopt the 
double-Gaussian function as an emission line profile model. 
Note that our results do not significantly depend on the choice 
of the model since the difference in flux measurements is ~ 
5%. 

In the rest-frame UV spectra, many AGN emission lines are 
blended; Lya A1216+N VA1240, Si IVA1397+0 iv]A1402, 
and He IIA1640+O III]A1663+A1 IIA1671. Thus, it is nec- 
essary to perform a multi-component fitting analysis for se- 
cure flux measurements. We simultaneously fitted all 10 
emission lines, that were used as BLR metallicity indica- 
tors. First, we divided these emission lines into two groups 
based on their ionization degree, and assumed that the emis- 
sion lines in each g roup have the same velocity profile (see 
INagao et aT1l2006bl) . Specifically, we categorized N VA1240, 
O iv]A1402, N iv]A1486, C IVA1549 and He HA1640) in 
the high-ionization group, while Si IIA1263, Si IVA1397, 
O Hi] A 1663 and Al I I A 1671 in the low-ionization group, fol- 
lowing |Nagaoetal] d2006bl) . We adopted the same velocity 
width and the velocity shift for each group. We excluded the 
spectral range between 1570A and 1631 A from the fit, since 
an u nidentified emission feature is reported in th i s range (see, 
e.g..lWilkeslll98l lBovIe1[T990l: ILaor et al.lfl994l: INagao et alJ 
I2006bl) . In the case of Lya, the absorption by the intergalac- 
tic matter (IGM) affects the line profile significantly, partic- 
ularly below 1210A. Thus, we treated Lya separately, by al- 
lowing the velocity dispersion and velocity shift to be free. 
We excluded the spectral range below 1210A from the fitting 
procedure, because of the IGM absorption. For the contin- 
uum fitting, we used a power-law and determined the slope 
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FIG. 2. — Comparison of different fitting functions for the C IV A1549 emis- 
sion line of PG 0003+158. Double-Gaussian (upper left), Gauss-Hermitian 
(upper right), modified-Lorentzian (lower left) and 2 power-law (lower right) 
functions are examined. The red lines denote the fitting results, and the blue 
line in the upper left panel represents each Gaussian component. Residual 
spectrum is shown at the bottom in each panel. Masked regions are indicated 
with gray hatches. 
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FIG. 3. — Examples of Multi-component fitting analysis. Each panel 
shows Lya A1216+N VA1240 (upper left), Si IVA1397+0 IV]A1402 
(upper right), N IV]A1486+C IVA1549 (bottom left), and 
He IIA1640+O III]A1663+A1 IIA1671 (bottom right), respectively. 
The color of lines are the same as in Figure 2. The dashed lines indicate the 
center of each emission line. 

by using three spectral windows (1345 A - 1355 A, 1445 A - 
1455A, and 1687A - 1697 A), where no strong emission lines 
are present. 

Figure 3 shows examples of multi-component fit- 
ting for 4 different spectral regions: from upper 
left to bottom right, 1) Lya A1216+N VA1240, 2) 
Si IVA1397+0 iv]A1402, 3) N iv]A1486+C IVA1549, 
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FIG. 4. — Fitting results of 2 representative objects: PG 0003+018 with a 
high quality COS spectrum (top) and PG 0007+106 with a low quality IUE 
spectrum (bottom). The inset panels show the fit for weak lines, namely, 
Si IV+O IV] (left) and He II+O III]+A1 II (right). No inset panel means that 
fluxes of these weak lines were not measured due to the low signal-to-noise 
ratio. The green dashed lines in the inset panels represent the center of each 
line as shown in Figure 3. The fitting residual is shown in the lower panel for 
each object. The color of lines and masked regions are the same as in Figure 
2. 

and 4) He HA1640+O iii]A1663+A1 nA1671. The line center 
is denoted with dashed lines. 

3.2. Fitting result 

Based on the multi-component fitting using double- 
Gaussian profiles, we measured the flux of 10 broad emission 
lines in the rest-frame UV spectra for the sample. In Fig- 
ure 4, we present the fitting results of 2 representative targets; 
PG 0003+199 with high quality COS data and PG 0007+106 
with low quality IUE data. In some cases, spectral quality 
is too low to fit the weak emission lines, i.e., N VA1240, 
Si IVA1397, O iv]A1402, C IVA1549, and He HA1640, thus 
we only measure the flux of N V and C IV. The weak lines 
(i.e., Si IV+O IV] and He II+O Ill]+Al II) are shown in the 
inset panels only when these lines were successfully fitted. In 
summary, we measured the N V and C IV fluxes for the en- 
tire sample (70 objects) while we measured the flux of Si IV, 
O IV], He II for a subsample of 34 objects. Table 2 lists the 
measured fluxes and the inferred uncertainties, which were 
estimated by averaging the signal-to-noise ratio of each pixel 
within the spectral range of each lines. We list the sum of 
Si IVA1397 and O IV]A1402 fluxes instead of individual flux 
measurements, since the sum of two lines will be used to com- 
pare with the flux of the C IV line. 

4. RESULT 
4.1. Comparison of emission-line fluxes 

In this study, we use three line flux ratios as metallicity indi- 
cators among various diagnostics previously proposed, due to 
the limited data quality; namely, N VA1240/C IVA1549, 
(Si IVA1397+0 iv]A1402)/C IVA1549, and 
N vA1240/He IIA1640 (hereafter N v/C IV, 
(Si IV+O IV])/C IV, and N v/He II, respectively). Be- 
fore using them as metallicity indicators in comparison with 
AGN properties, here we first examine the correlation among 
the emission line fluxes. 

Figure 5 (top panel) presents the comparison between the 
C IVA1549 and He HA1640 fluxes, which are used as a 
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TABLE 2 

Measurements of emission Line fluxes 



Object NV SilV+OIV] C TV Hell 

(10" 14 ergs"') 

(1) (2) (3) (4) (5) 



NOTE. — Col. (1): Target ID. Col. (2): line flux and error of N V. Col. (3): line flux 
and eiTor of Si IV+O IV. Col. (4): line flux and error of C IV. Col. (5): line flux and 
error of He II. 



denominator of metallicity diagnostics. As photoionization 
models predict that the flux ratio of th ese two lines does n ot 
depend on Zblr (see, e.g., Figure 29 in lNagao et al.j|2006bl) . a 
clear linear relation with a small scatter (0.20 dex) is present, 
showing that the flux ratio of these two lines is nearly con- 
stant. The constant flux ratio between C IV and He II suggests 
that our flux measurements of the weak He II line are reason- 
able although the He II is difficult to fit due to blending with 
Al II. 

In Figure 5 (bottom panel) we compare the N VA1240 flux 
and the sum of the Si IVA1397 and O IV]A1402 fluxes, which 
are used as a numerator of metallicity indicators. Again a 
clear positive correlation is present between them, with a 
somewhat larger scatter than that shown in the comparison 
between the C IV and He II fluxes. This larger scatter is partly 
caused by the fact that the N V flux depends on Zblr as well 
as the relative abundanc e of N V (see, e.g.. iMatsuoka et ail 
1201 lbtlAraki et al.l2012l) , while the sum of the Si IV and O IV] 
fluxes mainly depends on Zblr- 

4.2. Comparison among metallicity indicators 

We compare 3 line flux ratios, namely, N v/C IV, 
(Si IV+O IV])/C IV, and N v/He II, as metallicity indicators 
adopted in this study. As shown in Figure 6, the flux ratios 
of N v/C IV and N v/He II present a clear correlation with 
a relatively small scatter, reflecting the constant flux ratio be- 
tween C IV and He II (see Figure 5). The comparison between 
(Si IV+O IV])/C IV and N v/C IV also shows a relatively good 
correlation although a few outliers dominate the scatter. In 
the case of (Si iv+O iv])/C iv and N v/He n, the comparison 
shows a less clear correlation, probably due to the larger com- 
bined uncertainties on the flux measurements of weak lines 
(Si iv+O iv] and He n). 

4.3. Comparison between metallicity and AGN properties 

In this section, we investigate the correlation between BLR 
metallicity inferred from the emission line ratios and AGN 
properties, i.e., black hole mass, luminosity, and Eddington 
ratio, using the selected low-z QSOs. Figure 7 compares 3 
metallicity indicators (N v/C IV, (Si IV+O IV])/C IV, and 
N V/He II) with AGN properties (M B h, L boh and L hol /L EM ), 
respectively. A positive correlation is present between 
N V/C IV or N V/He II with AGN luminosity, indicating 
a luminosity dependence of the BLR metallicity, while it 
is less certain in the case of (Si IV+O IV])/C IV, probably 
due to the larger measurement uncertainty in Si IV+O IV] 
flux. The apparent luminosity - BLR metallicity correlation 
of low redshift QSOs is similar to the results from previous 
studies based on high redshift AGNs (IShemmer et al. l l2004t 
IWarner et al.ll2004llNagao etaill2006bl) . 

We investigate which parameter between Mbh and Lboi/^Edd 
is more fundamental in driving the observed Lagn _ Zblr re- 
lation. As for Mbh, there is no obvious Mbh dependence on 
N V/C IV and (Si IV+O IV])/C IV while there is a possible 
positive correlation between N v/He II and Mbh- These be- 
haviors are in contrast to previous results obtained at high red- 
shifts, where significant positive correlations were reported 
between the metallicity ind icators and Mbh (IWarner et al.l 
l2004tlMatsuoka et al.ll20TTbh . 

On the contrary, clear positive dependences on Lboi/^Edd are 
present in N v/C IV and N v/He II, while it is less certain in 
the case of (Si IV+O IV])/C IV, due to the lack of objects 
in the range of Lboi/^Edd < -1-5. For these low Eddington 
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PG 1202+281 


8.3 ±0.8 


— 


72.9 ±5.1 


— 


PG1211+143 


80.6 ±3.3 


27.5 ±1.0 


171.0± 12.1 


— 


PG12 16+069 


50.8 ±8.7 


19.3±2.1 


105.6±9.6 


10.6 ±0.9 


PG 1226+023 


168.0±17.1 


59.8 ±3.0 


305.5 ±13.5 


32.3 ±1.4 


PG 1229+204 


15.0±0.6 


47.5 ±2.0 


156.4±4.8 


11.4±0.3 


PG 1244+026 


3.7 ±0.3 


— 


11.4± 1.2 


— 


PG1259+593 


25.5 ±6.0 


13.9±2.0 


28.1 ±4.3 


2.4 ±0.3 


PG1302-102 


46.8 ±1.0 


— 


52.4 ±0.6 


— 


PG 1307+085 


64.3 ±4.2 


— 


114.6±4.6 


— 


PGI310-108 


36.4 ±3.4 


— 


140.3 ±15.2 


26.1 ±2.4 


PG1322+659 


31.1 ±1.2 


— 


51.4±2.1 


6.6 ±0.3 


PG1341+258 


22.4 ±2.6 


— 


34.7 ±5.3 


— 


PG1351+695 


18.1 ± 1.8 


10.9 ±0.7 


99.6± 11.1 


5.3±0.9 


PG1352+183 


29.9 ±2.6 


— 


47.4 ±2.8 


— 


PG1402+261 


15.0±1.0 


— 


68.9 ±3.4 


— 


PG1404+226 


9.5 ±1.3 


— 


13.1 ±1.6 


— 


PG1415+451 


44.0 ±2.7 


20.4±1.1 


55.0±2.8 


6.6 ±0.3 


PG1416-129 


6.7 ±0.9 


— 


66.4 ±4.3 


— 


PG1425+267 


8.6±0.4 


— 


49.4 ±1.5 


— 


PG1426+015 


94.3 ±6.8 


— 


260.7 ±36.8 


— 


PG1427+480 


13.8±0.7 


9.3 ±0.4 


43.1 ±2.0 


6.1 ±0.3 


PG1434+590 


186.9 ±4.2 


86.7 ± 1.4 


379.3 ±16.2 


— 


PG1435-067 


32.3 ±2.6 


— 


62.4 ±4.1 


— 


PG1440+356 


96.1 ±5.1 


54.7 ±2.2 


119.1 ±4.9 


20.9 ±0.7 


PG1444+407 


36.0±2.1 


— 


32.6 ±2.8 


— 


PG1448+273 


10.0 ±0.8 


3.7 ±0.2 


13.2± 1.5 


4.1 ±0.5 


PG1501+106 


18.0±0.5 


— 


219.0±5.6 


9.8 ±0.2 


PG15 12+370 


32.9 ±5.8 


8.3 ±1.0 


70.2 ±7.3 


4.5 ±0.8 


PG15 19+226 


22.4 ±3.4 


29.0 ±6.0 


49.0 ±5.9 




PG1534+580 


40.8 ± 1.5 


25.1 ±0.5 


187.5 ±6.0 


14.3 ±0.6 


PG1543+489 


25.7 ±0.9 


11.8±0.3 


23.7 ±0.5 




PG 1545+2 10 


36.9 ±2.0 




101.6±5.4 


6.9 ±0.4 


PG1552+085 


12.7 ±2.2 




30.7 ±4.4 




PG1612+261 


11.0± 1.3 




67.8±5.1 


2.5±0.1 


PG1613+658 


63.7 ± 1.2 


19.2 ±0.4 


175.7±4.5 




PG1617+175 


12.5 ± 1.2 




41.2±2.1 




PG1626+554 


55.0±3.0 


18.5 ± 1.0 


78.9 ±3.8 


5.7±0.3 


PG2 112+059 


23.8 ±3.8 


9.1 ± 1.1 


16.0±2.0 




PG2 130+099 


34.9 ±0.9 




120.8 ±6.7 


14.9 ±0.9 


PG2214+139 


9.8 ±2.7 




88.8 ±14.3 


8.0 ± 1.5 


PG2233+134 


10.4 ±0.8 


4.4 ±0.2 


8.0 ±0.3 




PG2251+113 


21.0±2.2 


5.3±0.2 


29.4 ±1.2 




PG2304+042 


12.5 ± 1.6 




89.2±8.8 




PG2308+098 


23.5 ±4.5 


6.4 ±0.8 


37.0±4.8 
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FIG. 5. — Comparison of emission-line fluxes. Top: Comparison between 
the He IIA1640 and C IVA1549 fluxes. Bottom: Comparison between the 
sum of the Si IVA1397 and O IV]A1402 fluxes and N VA1240 flux. The 
derived Spearman's rank order correlation coefficients and their statistical 
significance are +0.795, 3.2 X 10~ 8 (top) and +0.751, 4.5 X 10~ 7 (bottom). 



ratio objects, Si IV+O IV] lines are very weak and no secure 
measurements are available. 

In order to examine the statistical significance of these pos- 
sible correlations, we applied the Spearman's rank-order test 
to the data. The derived Spearman's rank-order correlation 
coefficients (rs) and their statistical significance (p), which 
is the probability of the data being consistent with the null 
hypothesis that the flux ratio is not correlated with an AGN 
parameter, are given in Table 3. The rank-order tests suggest 
that there are statistically significant positive correlations of 
N v/C IV with Lboi and Lboi/^Edd, and a positive correlation of 
N v/He II with Lboi- In contrast, all three diagnostic flux ratios 
show no statistically significant correlation with Mbh- 

5. DISCUSSION 

We compare our results obtained for AGNs at z < 0.5 
with the previous resu lts obtained for high-z QSOs (z-~2.5; 
iMatsuoka et al.ll2"01 lbl) . in order to investigate possible red- 
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FIG. 6. — The relation among three metallicity indicators. Red line rep- 
resents error weighted linear fit to the data. The slopes of the linear fit and 
the data dispersion are shown at the upper-left comer in each panel. The 
derived Spearman's rank order correlation coefficients and their statistical 
significance are +0.726, 2.0 X lO" 6 (top), +0.601, 2.1 X 10~ 4 (middle) and 
+0.329, 1.6 X 10-' (bottom). 
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FIG. 7.— 
or mass bin 



6 7 8 9 10 

Log (M BH /M Q ) 

The relation between metallicity indicators 
as indicated in the top panels. 



44 45 46 47 -3 -2 -1 

Log L bol (erg/s) Log (L/L Ed<j ) 

and AGN properties (A/bh, £bol, and Lbol/^Edd)- The symbols and colors represent the Eddington ratio 



shift evolution and the origin of the metallicity scaling rela- 
tions in AGNs. In Figure 8, we overplot the metallicity indica- 
tors a s functions of Mbh and L /L^ for high-z QSOs adopted 
from dMatsuoka et al.ll201 lbl) . along with the measurements 
of low-z QSOs. The low-z objects are more dispersed than 
high-z objects, partly because of the larger measurement un- 
certainties of the low-z objects. Note that the emission line 
fluxes were measured for individual objects in the low-z sam- 



ple, while the composite spectra of high-z QSOs in each black 
hole mass and Eddington ra tio bin were used for the emission 
line flux measurements (see Mats uoka et al.ll201 lbl) . 

Interestingly, the low-z and high-z samples appear to show 
different trends in Figure 8 (top panels). As mentioned in 
§4.3, the low-z sample shows only marginal correlations be- 
tween metallicity indicators and Mbh, while the high-z sam- 
ple clearly shows positive correlations of metallicity indica- 
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TABLE 3 

Results of Spearman' s rank-order correlation test 



Flux ratio 




Mbh 




Lbol 


ibol/^Edd 


(1) 




(2) 




(3) 


(4) 


N V/C IV 


'"s 


= +0.076 




= +0.487 


r s = +0.473 




P = 


5.3 x 10-' 


/' = 


2.1 x 10~ 5 


p = 3.6 x 10~ 5 


(Si IV+O IV])/C IV 


rs 


= -0.128 


rs 


= +0.038 


r s =+0.381 




P = 


4.8 X 10-' 


P = 


8.4 x 10-' 


p = 2.9 X 10~ 2 


N v/He II 


rs 


= +0.417 


rs 


= +0.622 


r s =+0.311 




P = 


1.6 x 10~ 2 


P = 


1.1 x 10~ 4 


p = 7.8 X 10~ 2 



* -0.4*log(L/L Eoa ) 

* -0.8^log(L/L E0a )<-0.4 

* -1.2^log(L/L Efla )<-0.8 
■ -1.6^log(L/L Efla )<-1.2 

* log(L/L Efla )<-1.6 



★ 9.4^log(M ah /U Q ) 

• 8.8glog(M ah /M )<9.4 

* 8.2glog(M aH /U )<8.8 
■ 7.6glog(M aH /U )<8.2 
X log(U, H /U„)<7.6 




Log (M BH /M S ) 



-2 -1 

Log (L/L E00 ) 



FIG. 8. — Top: Metallicity indicator (N V/C IV flux ratio) as functions 
of Mbh (left) and L/LEdd (right). Low-z QSOs are repres ented with vari- 
ous sy mbols while high-z QSOs are denoted by solid lines iMatsuoka et al.l 
1201 lbl) . The colors represent different mass and Eddington ratios as indicated 
in the upper panels. Bottom: Metallicity indicator vs. Mbh (left) and L/Le<m 
(right), after excluding low Mbh and low Eddington AGNs from low-z sam- 
ple for proper comparison. 

tors with Mbh- On the other hand, the metallicity indicators 
of the low-z sample show much stronger positive correlations 
with L/Lndd than with Mbh, while the correlations between 
metallicity indicators and L/Lsdd are less evident in the high- 
Z sample. 

The difference in the trend with Eddington ratios between 
low-z and high-z samples is partly caused by the much wider 
range of L/Lem ( _ 3 < log[L/LEdd] < 0) covered by the low- 
z sample than that covered by the high-z sample (-1 < 
log [i/Lp^d] < 0). As shown in Figure 8, low-z AGNs with 
log[L/LEdd] < -1-5 (that is not covered in the high-z sample) 
show systematically lower N v/C IV flux ratios, leading to 
the more evident dependence on L/LEdd- Therefore, our re- 
sults do not necessarily suggest that the L/L E dd dependence 
on the metallicity indicators is systematically different be- 
tween low-z and high-z QSOs. For proper comparison, low 
Eddington AGNs (log[L/LEdd] < -1 -5) are required at high-z. 
These results imply that the accretion activity of black holes 
are closely related with metal enrichment at the central part of 
host galaxies. 



The dependence of metallicity indicators on Mbh appears to 
be different between the low-z and high-z samples. At high- 
z, more massive AGNs have higher metallicity although the 
Mbh range is small (8.5 < log Mbh< 10). In contrast, the 
low-z sample shows much larger scatter without clear trend 
between Mbh and metallicity indicators. However, the large 
scatter in the low-z sample is partly caused by the low Edding- 
ton AGNs, which have systematically lower N v/C IV flux ra- 
tios, as described as the metallicity-Eddington ratio relation. 
However, if AGNs with similar Eddington ratios are selected 
(e.g., same color objects in Figure 8), then there appears to 
be a weak trend of metallicity with Mbh, suggesting that at 
low-z, Zblr has weak dependency on Mbh at fixed Eddington 
ratios. 

For high-z QSOs, the tightness of the observed Zblr -Mbh 
correlation is probably caused by the limited range of the Ed- 
dington ratios since only high Eddington QSOs were included 
in the high-z sample. However, inclusion of low Edding- 
ton AGNs will presumably weaken the correlation as in the 
case of the low-z sample. As a consistency check, we match 
the ranges of Mbh and Eddington ratio between high-z and 
low-z samples, by excluding low Mbh and low Eddington ob- 
jects from the low-z sample as presented in the bottom panels 
of Figure 8. As expected, both high-z and low-z samples in 
the matched dynamical range show similar metallicity depen- 
dency on both Mbh and Eddington ratio although the Mbh- 
Zblr relation is much weaker in low-z than at in high-z. 

Our results imply that the BLR metallicity of low-z AGNs 
mainly depends on the Eddington ratio, and weakly depends 
on Mbh- Currently, it is unknown how Zblr scales with av- 
erage gas metallicity of host galaxies. Nevertheless, the ob- 
served metallicity dependency on Mbh may imply that there 
is connection among BH growth, gas enrichment, and galaxy 
evolution. 

Assuming that Z B lr correlates with gas metallicity of host 
galaxies, we try to understand the observed relations with sev- 
eral scenarios. The correlation between Mbh and Zblr at 
high-z can be interpreted as a consequence of the combina- 
tion of the gala xy mass-metallicity relation and the Mbh ~ 
Mbuige relation (IWarner et al.l 120031: iMatsuoka et al.l l2011bb . 
For QSOs at z ~ 2 - 3 that corresponds to the peak of the 
quasar activity in the cosmological timescale, it has been of- 
ten claimed that the major merger triggers the AG N activity 
through the efficient mass fueling onto black holes dHasingerl 
120081: iLi et al.ll20Tob . Here the major merger event reduces 
the angular momentum of gas clouds at the nuclear regions 
of the quasar host galaxies, providing efficient mass fueling. 
In this case, the metallicity of accreting gas onto the nucleus 
may be characterized by the mass-metallicity relation of the 
host galaxy. Regarding the black hole to host galaxy connec- 
tion, the Mbh -Mbuige relation has not been observationally 
defined at high-z although several studies in dicated that Mbh- 
to-Mb u i e e ratio increases with redshift (e.q.. | Woo et al.1 120061 
120081: iBennert et al1l20H l20llt iMerloni et al.ll2010l) . Thus, 
if we assume a scaling relation between Mbh and galaxy mass, 
presumably with higher normalization than the local Mbh- 
Mbuige relation, then black hole mass scales with both stellar 
mass and gas metallicity, hence the Mbh-Zblr relation is nat- 
urally expected at high-z. 

The weaker correlation of Mbh with Zblr at low-z can be 
interpreted as combination of two additional effects. First, gas 
metallicity increased from high-z to low-z by further star for- 
mation after major BH growth. In other words, for the same 
galaxy mass (or black hole mass), metallicity has been in- 
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creased, particularly at the nuclear region, leading to higher 
metallicity at fixed Mbh compared to high-z objects. For 
example, nuclear star formation induced by secular process 
(e.g., bar instability) and galaxy interaction m ay sufficiently 
increase metallicity (e.q., |Maiolino et aT1l2008l) . Another dif- 
ference between low-z and high-z is the gas fraction, espe- 
cially for massive galaxies like quasar host galaxies. The 
minor merger process is mor e important as an AGN trigger- 
ing m echanism at low-z (e.q.. lTaniguchill999l:ICisternas et al] 
l201lh . Thus, as a consequence of minor mergers between 
a gas-poor massive galaxy (i.e., a quasar host galaxy) and a 
relatively gas-rich less-massive galaxy, the metallicity of ac- 
creting gas onto the nucleus can be largely affected by the 
chemical property of the merging, less-massive galaxy. In this 
case, the metallicity of accreting gas does not simply reflect 
the mass-metallicity relation of the host galaxy. Second, the 
MBH-to-Mbuige ratio may change over cosmic time. If black 
hole mass was higher at fixed stellar mass at high-z, then by 
combining galaxy mass-metallicity relation and MBH-Mb u i ge 
relation, we may expect similar MBH-metallicity relation in 
low-z, but with a different normalization. In other words, at 
fixed metallicity, Mbh is lower in low-z than in high-z. Al- 
though it is beyond the scope of the current work to quantify 
and compare these two effects, it is reasonable to conclude 
that coupling between black hole mass with gas metallicity 
becomes much weaker with decreasing redshift. 

6. SUMMARY & CONCLUSION 

To investigate the chemical properties of low-z QSOs, we 
measured the flux ratios of the rest-frame UV emission lines 



as metallicity indicators using a sample of 70 low-z PG QSOs 
at z < 0.5. By comparing BLR gas metallicity with black hole 
mass, luminosity and Eddington ratio, we find that Zblr cor- 
relates with Eddington ratio while Zblr shows much weaker 
correlation with Mbh, indicating that the metal enrichment at 
the central part of host galaxies is closely connected to the ac- 
cretion activity of AGN. These trends are different from high- 
z QSOs, which shows a tighter correlation between Zblr and 
Mbh and a weaker correlation between Zblr and Eddington 
ratio. The apparent difference between low-z and high-z sam- 
ples seems to be caused by the limited dynamical range in 
the high-z sample. Various star formation mechanism can in- 
crease BLR gas metallicity the cosmic time, increasing the 
scatter in the metallicity correlation with properties of AGN 
in low-z. 
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